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Abstract 

I show that the so-called causality paradox of time-dependent density functional theory arises 
from an incorrect formulation of the variational principle for the time evolution of the density. The 
correct formulation not only resolves the paradox in real time, but also leads to a new expression 
for the causal exchange-correlation kernel in terms of Berry curvature. Furthermore, I show that 
all the results that were previously derived from symmetries of the action functional remain valid in 
the present formulation. Finally, I develop a model functional theory which explicitly demonstrates 
the workings of the new formulation. 
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I. INTRODUCTION 

Time-dependent density functional theory (TDDFT) is becoming a standard 

tool for the computation of time-dependent phenomena in condensed matter physics and 
quantum chemistry. Naturally the growing number of application has generated a new 
interest in the foundations of the theory (see, for example the recent critique by Schirmer and 
Drew [4], and the rebuttal by Maitra, Burke, and van Leeuwen jlj]). In this paper I address 
the so-called "causality paradox" , a problem that has troubled TDDFT for many years [6] , 
and has been the object of many discussions and technically sophisticated resolutions Q,u3, 

[9I, [hJ. I do not disagree with those resolutions, but I wish to propose a new one, which 
I find technically simpler, more direct, and closer to the spirit of the original formulation of 
TDDFT. 

As an introduction to theproblem, let us recall that the formal basis of TDDFT is the 
Runge-Gross (RG) theorem[2j, which establishes a biunivocal correspondence between the 
time dependent particle density n(r, t) of a many-body system and the potential v(r, t) 
that gives rise to that density starting from an assigned quantum state |^o) at the initial 
time t = 0. According to the RG theorem the potential that gives rise to n(r, t) starting 
from j^o) is determined by n(r,t) and |^o) U P to a time-dependent constant. Similarly, 
the time-dependent quantum state \ip(t)) is determined by n(r,t) and j^o) U P to a time- 
dependent phase factor. In this sense, both v(r, t) and \ijj{t)) are Junctionals of n(r, t) and 
|^o) over a time interval < t < T. They should be denoted by v[n, |V>o) j r, t] and \ip[n, \ipo)}) 
respectively. From now on, however, the dependence on the initial state will not be explicitly 
noted and we will write simply v[n; r, t] and | *0 [ri] ) . Admittedly, there is no proof that every 
reasonable density can be produced by some local potential: but it is generally assumed 
that the densities for which the theorem holds are dense enough in the space of densities to 
provide an arbitrarily good approximation to the physical densities one might encounter in 
real life. In this paper I will assume tout court that all the time-dependent densities can be 
produced by some local potential: i.e., all time-dependent densities are v-representable. 

Another theorem, proved by van Leeuwen in 1999, ll| extends and strengthens the RG 
theorem. According to van Leeuwen's theorem, the density n(r, t) which evolves in an 
interacting system under the action of an external potential v(r,t) starting from an initial 
state \ipo), can be reproduced in a noninteracting system evolving under the action of an 
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appropriate and uniquely determined potential v s (r,t), starting from any initial state \ip s o) 
that has the same density and divergence of the current density as \ipo). This theorem 
provides the basis for the extremely useful Kohn-Sham method of calculating the density. 
The effective potential v s (r,t) - a functional of n(r, t), \ipo), and \ip s o) - is known as the 
Kohn-Sham potential. The difference v(r, t)—v s (r, t) — vu(r, t), where vh{t, t) is the Hartree 
potential, is known as the exchange-correlation (xc) potential, denoted by v xc [n;r,t] - also 
a functional of n(r,t), \ipo), and l^o)- 

In general the potentials v(r, t), v s (r,t) and v xc (r,t) depend on the density n(r',t') at 
different positions and earlier times t! < t, but cannot be affected by changes in the density at 
later times t! > t. This obvious causality requirement implies that the functional derivatives 
of these potentials with respect to n(r,t') and, in particular, the exchange-correlation kernel 
/ xc (r, t; r', t') = 5v xc [n] r, t]/8n(r', t') vanish for t < t'. 

An interesting question is whether the potentials v(r,t), v s (r,t) and v xc (r,t) can be 
generated from functional derivatives of an action functional A[n, \ipo)} (denoted from now 
on simply as A[n]) with respect to the density, in close analogy with static DFT, where 
the potentials are functional derivatives of energy functionals with respect to the density. 
The existence of such a representation was sug gested by RG in their original paper [2fl, and 



was subsequently used by this author 12|, [l3| to derive several theorems in TDFT. In the 



mid-nineties, however, it became clear that the representation was problematic to say the 
least [6]. If the potential could be written as a functional derivative of an action functional, 

v\n-vt]= 6A[n] (1) 

then we should also have 

5v[n;r, t] 5 2 A[n] 



(2) 

6n(r',t') 8n{r,t)8n(r',t') v ; 

But this equation is patently false, because the left hand side is different from zero only 

for t > t' (by the causality requirement), while the right hand side is symmetric under 

interchange of t and t' . 

This startling observation became quickly known as the causality paradox and prompted 

re best known is the Van Leeuwen's con- 

n 



several sophisticated resolutions 
struction of a "Keldysh action" in pseudotime [s 



9(. More recently Mukamel 10| has shown 



how to construct causal response functions from symmetrical functional derivatives corre- 
sponding to "Liouville space pathways" . The gist of these resolutions is that causality is not 



violated, but one must use a more abstract mathematical apparatus (Keldysh formalism, or 
the Liouville superoperator method) in order to connect functional derivatives of the action 
to causal response functions. 

In this paper I re-examine the "paradox" from a more elementary point of view. I show 
that the variational principle for the time-evolution of the wave function, when properly 
implemented as a variational principle for the density, yields an expression for the potential 
as the sum of two terms: (1) the functional derivative of the RG action and (2) a correction 
term, which cannot be expressed as a functional derivative, but is still simple enough to 
be included in all the formal proofs. So the gist of the present resolution is that we learn 
to write the potential as functional derivative of an action plus a boundary term. Among 
other benefits, this approach explains why theorems that were originally proved under the 
incorrect assumption ([TJ) [l2, Q, turned out to be true after all. Furthermore, it leads to 
interesting expressions for the inverse of the density-density response function and the xc 
kernel in terms of "Berry curvature" . 

This paper is organized as follows. In the next section I discuss in detail the failure of 
the stationary action principle for the density, and show how the correct causal expressions 
for the xc potential and the xc kernel are derived from a modified variational principle. 
In section II III I explain why in many cases one can still pretend that the xc potential is 
the functional derivative of the xc action and get correct results. Appendix |A] clears up a 
technical point about the equal-time singularities of causal response functions. Finally, Ap- 
pendix [B] presents a pedagogical "time-dependent position density functional theory", which 
is conceptually equivalent to the full-fledged TDDFT but can be solved exactly, illustrating 
the workings of the new formulation. 



II. VARIATIONAL PRINCIPLE FOR THE DENSITY 

The starting point is the time-dependent quantum variational principle, [y| according to 
which the time- dependent Schrodinger equation is equivalent to the requirement that the 
action 

M\i>)]= f (i/>(t)\idt-H v \i>(t))dt (3) 
Jo 
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(h = 1) be stationary {5 Ay = 0) with respect to a arbitrary variations of the wave function 
which vanish at the ends of the time interval < t < T, i.e. \5ip(0)) = \5ip(T)) = 0. Here 



H V = H + j V(r,t)n(r)dT (4) 

is the sum of the internal hamiltonian H (kinetic + potential) and the interaction with an 
external time-dependent potential field V(r,t). 

The proof is straightforward. The variation of Ay induced by a variation 5\ip) = \5ip) is 

5A V U)] = f {5^mdt-Hv\m)dt 
Jo 

+ / (i>(t)\id t -Hv\Sm)dt , (5) 
Jo 

and the second term on the right hand side can be integrated by parts to yield 

6A v [\il>)] = f \srl>(t)\id t -H v \il>(t))dt 
Jo 

+ / ((id t -H v )m\m))dt 

Jo 

+ i(m\m)) to ■ (6) 

The last term on the right hand side vanishes by virtue of the boundary conditions on \Sip(t)), 
and the vanishing of the first two terms is equivalent to the time-dependent Schrodinger 
equation (id t — Hv)\tp(t)) = 0. 

Since \ip(t)) is, by virtue of the RG theorem, a functional of n and IV'o)) Runge and Gross 
suggested that a stationary action principle for the density could be formulated in terms of 
the functional 

rp 

A v [n] = j {rj)[n\\idt- H\rj)[n\)dt 
Jo 

= A [n\ - / V(r,t)n(r,t)drdt , 
Jo 

(7) 

where the "internal action" 

A [n} = [ {il>[n]\idt - H \il;[n])dt (8) 
Jo 

is a universal functional of the density and the initial state. Then, setting 5A V = for 
arbitrary variations of the density we easily find 



when n(r, t) is the density corresponding to V(r, t). This implies that the external potential, 
viewed as a functional of the density, is the functional derivative of the internal action with 
respect to the density: 

^•^iMr (10) 

and the time evolution of the density is determined by requiring v[n;r,t] = V(r, t), where 
V(r, t) is the actual external potential. The only problem with Eq. (jlOj) . which would other- 
wise be very useful, is that it plainly contradicts causality, as discussed in the introduction. 
What went wrong? 

The problem arises from the fact that the Frenkel variational principle 5 Ay = is valid 
only for variations of that vanish at the endpoints of the time interval under considera- 
tion, i.e. at t — and t = T. But a variation of the density at any time t < T inevitably 
causes a change in the quantum state at time T. Therefore we can only set \5ip(0)) = 0, but 
have no right to set \Sip(T)) = 0. Taking this into account, and going back to Eq. ([6]) we see 
that the correct formulation of the variational principle for the density is not 5 Ay = but 

M v [n]=tWn]|^[n]) . (11) 

Here I^tH) = l^fajT]) * s the quantum state at time T regarded as a functional of the 
density (and of course of the initial state). So we see that the action functional is not 
stationary, but its variation must be equal to another functional of the density, which is 
given on the right hand side of Eq. ffTTT) . 

Taking the functional derivative of Eq. (iTTj) with respect to n(r, t) and making use of 
Eqs. (ED and ^TDJ) we get 

5i/j T [n] 



r n <L4 [n] I 



5n(r, t) 



(12) 



where 



j s a compact representation for the functional derivative of IV'tW) with 



5n(r,t) 

respect to density. 

This is the main result of this paper, since it shows that the external potential (and hence 
also the Kohn-Sham potential and the xc potential) is not merely a functional derivative 
of the Runge-Gross action Ao[n]. Notice that the additional "boundary term" is real, in 
spite of the i, because the quantum state IV'tM) is normalized to 1 independent of density, 
implying that 



ip T [n] 



5ip T [n] \ / 6if} T [n] 

5n(r, t) I \ 5n(r, t) 



Mn] ) (13) 



is a purely imaginary quantity. 

At first sight, however, Eq. (j!2p is still problematic because it appears to depend on the 
arbitrary upper limit of the time interval (T) and therefore also on the density at times 
if > t. However, this is only appearance. The point is that both the functional derivative of 
A [n] and the boundary term, considered separately, have a non-causal dependence on the 
density, but the dependence on n(t') with t' > t cancels out exactly when the two terms are 
combined! 

Let us show this in detail. Consider, for example, an increment of the upper limit of the 
time interval from T to T + AT. The density n(r, t) must be smoothly continued to the 
larger time interval [0, T + AT], and the quantum state satisfies in this time interval 

the Schrodinger equation 



{id t -H -v[n;t\Mn;t}) = , (14) 



where 



v[n;t) = J v[n\ r', t]n(r')dr' , (15) 

n(r) is the density operator, and v[n; r,t] is the potential that yields n(r,t). Now in view of 
Eq. (114jl the internal action over the time-domain [0,T] can be written as 

A [n] = r^tn;*']^^;*'])^^;*'])^' , (16) 



and its variation, due to the extension of the upper limit from T to T + AT is, to first order 
in AT given by 

AA [n] = AT(if} T [n]\v[n; T]\ip T [n\) . (17) 

Taking the functional derivative with respect to n(r,t), with t within the interval [0,T], we 
see that 

SMn] Am ,.,. 6v[n;T] ^ (lg) 



6v[n; T] 



Sn(r, t) 

The reason why we could take the functional derivative inside the expectation value on 
the right hand side of this equation is that (■0y[n]|u[n;T]|'02'[7i]) = f v[n; r', T]n(r', T)dr' 
depends on n(r,t) only through the potential functional v[n;r',T}: the density n(r',T) is, 



by definition, unaffected by a variation of the density at the earlier time t. [21] 

Consider now the variation of the boundary term of Eq. (fT2l) again due to the change of 
the upper limit of the time interval from T to T + AT. To first order in AT we have 

IV't+atN) - |VtN> = -iATH[n;t]\ip T [n]) , (19) 



where H[n, t] = Hq + v[n; T] is the full time- dependent hamiltonian regarded as a functional 
of the density. Substituting this in the variation of the boundary term we get 



—iA (ifa[n] 
ATU T [n) 



5n(r, t) 
H[n-T] 
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5n(r,t) 5n(r,t) 



H[n;T) 



i/) T [n] 



(20) 



This seemingly complicated expression contains the commutator between the Hamiltonian 
H[n;T] and the functional derivative with respect to n(r, t). This commutator is simply 
— ; which is evidently equal to — ^("'T) • Thus the variation of the boundary term is 
given by 



iA ( iPtW] 



5ip T [n] 



AT ( ipT [n] 



5v[n; T] 



Sn(r, t) 



(21) 



5n(r, t) i 

Combining the two equations (Fl8l) and fl2ip we see that net variation of the potential v[n; r, t] 
is exactly zero. This means that we have the freedom to change at will the upper limit T 
of the time interval: the value of v[n; r, t] will not change, as long as T remains larger than 
or equal to t. But this means that we can always choose T = t, and this proves that the 
potential at time t does not depend on what happens to the density at times later than t. 
QED. 

Now that we have the correct expression for the potential as a functional of the density 
it is easy to construct the other two potentials of interest, namely the Kohn-Sham potential 

and the xc potential. For the Kohn-Sham potential we simply have 

\8i> sT [n} 



5A 0s [n] , 
v s [n; r, t\ = - % { ipjr[n\ 



Sn(r, t) 



(22) 



Sn(r, t) 

where A 0s is the internal action for a noninteracting system and iV'arM) is the non 
interacting version of IV'tM) (starting from an initial state \ipso))- F° r the exchange- 
correlation potential we get 

SA xc [n] 



v xc [n; r,t] ■■ 
-i(ip T [n] 



5n(r, t) 



Sn(r, t) 



+ % ( ip sT [n) 



5i/; sT [n\ 



Sn(r, t) 



(23) 



where A xc [n] is the usual xc action defined in the Runge-Gross paper as the difference 
A [n] — A 0s [n] — where the last term is the Hartree action. Notice that v xc depends 

not only on the density, but also on the two initial states \ip ) and |-0 sO )- 
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The above formulas allow us to obtain an elegant expression for the functional derivative 
of the potential with respect to the density. Consider, for example, the functional derivative 
of v(n; r, t). From Eq. (TT21 we see that this is given by 



5v[n;r,t] 5 2 A [n] 
5n(r', f) ~ 8n(r,t)5n(r',t') 



ip T [n] 



8 2 ip T \n] 



— i 



Sip T [n] 



<Jn(r,t). 



.(24) 



6n(r, t)5n(r', if) / "\ 8n(r', V) 

The first two terms on the right hand side are symmetric under interchange of r, t and 
r',t', while the last term, which has the structure of a Berry curvature, is antisymmetric 
under the same interchange. Let us then subtract from Eq. (1241) the same equation with r, t 



and r', if interchanged. When t > if then Sv g™j* ^ ■ vanishes because of causality, so on the 



left hand side only 



8v [n;r,t] 



is left. And on the right hand side the symmetric terms cancel 



Sn{r',V) 

out, leaving only the Berry curvature term. The result is 



5v[n; r, t] 
5n(r',t') 



23m 



5n(r',t') 
5n(r',t') 



8ip T [n) 



5n(r, t) 

5ipT [n 



5ipr[n] 
Sn(r, t) 

(t > t'). 



Sipr[n] 



5n(r',t') 



(25) 



8n(r, t) 

Notice that the right hand side of this formula is independent of T (as long as T is larger 
that t and if) since, as we have shown, v[n; r,t] satisfies the causality requirements. 

The above argument determines f^fer^ for t>if, but leaves open the possibility of 
singular contributions at t — if . In Appendix [A] I show that indeed the functional derivative 
Sv[n; r, t]/5n(r', t') contains an equal-time singularity of the form 



C 6(t - t') + C x 8{t - t') + C 2 8(t - t') 



(26) 



where Co, C*i, and C 2 are functionals of the density at time t (= n t ) and functions of r and 
r'. 5 and 5 denote, respectively, the first and the second derivative of the 5-function with 
respect to its own argument. Thus the complete expression for Sv/Sn has the following form: 



Sv[n; r, t] 
8n(r', t>) 



26 (t - 1f)$tm 



6ipr[n] 
5n(r',t') 



5ip T [n] 



8n(r, t) 



+ S 00 [n t ;r,r')S(t-if) 



(27) 



where S'oof^t; r ; r '] is the differential operator C + Ci^ + C 2 ^. The coefficient of the leading 
term, C 2 , is independent of interactions, and the coefficient of the linear term, C±, vanishes 
in the linear response limit. An explicit demonstration of the equal-time singularities is 
provided in Eq. flB18j) of Appendix [B] 



Equal-time singularities also enter the expression of the exchange-correlation kernel. Tak- 
ing into account the fact that Ci is independent of interactions we find 

Sv xc [n; r,t] 



f xc [n;r,t,r',t'} = 



dn(r',t') 



5n(r't' 



5n(r, t) 



5n(r, t) I \ 5n(r', t') 

+ ACi[n t ; r, r')6(t - t') + f XC)00 [n t ] r, r')6(t - t'), 

(28) 

where ACi is the difference between the coefficients C\ in the interacting and non-interacting 
systems and f xc>0 o denotes the difference between the coefficients C in the interacting and 
non-interacting system. In the linear response regime, i.e. when the time- dependent po- 
tential is a weak perturbation to the ground-state, AC\ vanishes and f XC)CX > reduces to the 



well-known infinite-frequency xc kernel of linear response theory. [15| This behavior is demon- 
strated in Appendix [B] for our exactly solved model - see Eq. (IB21I) . 

Finally, I note that the adiabatic approximation to the xc kernel is given by 

ff c (r, t, r', ~ / SC)0 (r, r', t)5(t - if) , (29) 

where 



fxc,o{r, r, *) = f xc ,oo[n t ; r, r'] + / f xc [n; r, t, r', t']dt' 

Jo 



(30) 



is the integral of the exchange-correlation kernel over all times t' earlier than t. The implicit 
assumption here is that the retardation range of the xc kernel is much shorter than the time 
scale of variation of the density, so that the xc kernel can effectively be approximated as 
a 5-function on that time scale. The first term on the right hand side of this expression is 
the contribution of the "true" 5-function terms of Eq. (!28l) . The time integral in the second 
term is restricted to times strictly less than t. 

III. WHY DID "v xc [n;r,t\ = WORK? 

The incorrect representation of the xc potential as a functional derivative of the xc action 
played a significant role in the early development of TDDFT, particularly in the proof of 
theorems that depend on symmetries of the action functional. Consider, for example, the 
"zero-force theorem" [l2, 13], according to which the net force exerted by the xc potential 
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on the system is zero. This theorem was originally derived from the apparent invariance of 
the xc action under a homogeneous time-dependent translation of the density: 

A xc [n'\ = A xc [n] (31) 

where n'(r,t) = n(r + x(£),i), and x(i) is an arbitrary time-dependent displacement that 
vanishes at t — 0. The invariance of the action under this transformation implies 

6Axc[n] V r n(r)dv = , (32) 



5n(r, t) 

and an integration by parts leads to 

n^/A^dr = o . (33) 
v ; 6n(r,t) v ; 

This would be the zero-force theorem if we could identify g^(r ffi whh v xc [n; r,t] which, of 
course, is incorrect. 

Fortunately, the resolution of the puzzle is now at hand. The point is that we are making 
two errors, which are luckily compensating each other, leaving us with the correct result. 
The first error is in Eq. fl3Tl) : it is not true that the xc action remains invariant under the 



transformation n — > n'. The invariance of A xc was "derived" in Ref. (|12I|) by showing that 
the change of the internal action under this transformation depends only the density, not on 
the wave function, and therefore cancels out in the difference Aq — Aq s . However, we failed 
to include the boundary term i(ip T [n}\SipT[n]) , which does depend on the wave function and 
therefore does not cancel out, causing the xc action to vary, to first order in n' — n, by 

5 A xc = i (i) T [n] 1 5i) T [n\)-i (^t [n] \ 5^ sT [n] ) , (34) 



where = j^rfn']) — IV'tM) and I^stM) = IV'sT^']) — l^sT^])- Therefore, Eq. ( |3T|) 

must be amended as follows: 

A xc [n'] = A xc [n] + i(ip T [n}\Sip T [n]) - i(^ sT [n]\5ip sT , [n]) . (35) 

and Eq. ( l32|) is replaced by: 

5ip sT [n) 



SA xc [n] . , 



' x + i ( i> sT [n] 



5n(r, t) 



5n(r, t) 



V r n(r)rfr = . (36) 



Integrating by parts, and using the correct formula for v xc , Eq. (|23|) we do indeed recover 
the zero force theorem J n(r)V ' r v xc [n; r,t]dr = 0. 
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The lesson is quite general: we are allowed to pretend that the xc potential is the func- 
tional derivative of the action, provided we calculate that functional derivative incorrectly, 
i.e. ignoring the boundary contribution. This is exactly what we did (unwittingly) in our 
earlier papers. 

IV. CONCLUSION 

I believe that the foregoing analysis provides a straightforward and pedagogically trans- 
parent resolution of the causality paradox in TDDFT. Compared to the resolutions proposed 



in Refs. [8|, [10| the present approach is obviously much closer to the spirit of the original 
RG paper. Furthermore, our approach allows us to understand why in many cases we can 
get correct results from an incorrect representation of the xc potential. 

In closing I wish to emphasize that what we have derived here is a variational principle for 
the time-dependent density. A variational principle is not as strong as a minimum principle, 
yet it is strong enough to formulate a dynamical theory. While the absolute numerical value 
of the RG action has no physical meaning (because a multiplication of the wave function by 
an arbitrary phase factor changes its value by an arbitrary constant), it must be borne in 
mind that the action determines the dynamics through its variations, and those variations 
are independent of the arbitrary additive constant (a similar situation occurs in classical 
mechanics, since the Lagrangian is defined up to an arbitrary total derivative with respect 
to time). 
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APPENDIX A: EQUAL TIME SINGULARITIES IN 5v[n; r, t]/6n(r', t') 

Following Tokatly [l?j, Q] we write the exact local conservation laws for particle number 
and momentum: 

,9 t n + V-j = (Al) 

and 

mdtji + dj {mnUiUj + P^-) + ndiV = 0, (A2) 

where j is the current density, ji is its 2-th cartesian component, Ui = ji/n is the velocity field, 
dj denotes the derivative with respect to rj (with implied summation over repeated indices), 
and finally Py is the stress tensor. These equations are valid both for interacting and non- 
interacting systems and together define the time-dependent potential v as a functional of 
the density, provided the velocity field and the stress tensor are regarded as functionals of 
the density. Interaction effects enter implicitly throughf the form of these functionals. 

Taking the divergence of the second equation and making use of the first, we recast the 
system in the more explicit form 

di(ndiv) = md^n + mdi{uid t n) — m<9j(j ■ Viij) — didjPij. (A3) 

This equation can be formally solved, yielding 

v = G \md\n + mdi(uid t n) — m<9j(j ■ V-Uj) — didjPij] , (A4) 

where G is the inverse of the operator diudi. 

For the limited purpose of identifying the equal-time singularities in 5v/Sn we can ignore 
any retardation in the functional dependence of P^ and on the density. Then the right 
hand side of Eq. ( 1A4I) depends on the density and its first two derivatives d t n and d^n 
at time t. No higher derivatives are involved. Then taking the functional derivative with 
respect to n(r',t') we get a singularity proportional to S(t — t') from d^n(r,t), a singularity 
proportional to 5(t — t') from dtn(r, t), and, of course, a singularity proportional to 5(t — t') 
from the terms that do not contain time-derivatives of the density. 



y 



ree of interaction 



We can furthermore say that the coefficient of 5(t — t') is complete 
effects, since the interactions enter only in the functional P y - and U{. 22] And we observe 
that the S singularity vanishes in the linear response regime (small perturbations around the 
ground-state) because the current vanishes in the ground-state. 
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Very little can be said in general about the explicit form of the <5-function singularity. In 
the linear response regime, the dependence of P^ on density has been extensively studied, but 
only in local or semi-local approximations. 19J A fully nonlinear, but still local approximation 
to Py, known as nonlinear elastic local deformation approximation, has been formulated by 
Tokatly |l7| and studied by Ullrich and Tokatly 20] in a model calculation. An accessible 
review of this theory can be found in Chapter 8 of Ref . [l| • 



APPENDIX B: TIME-DEPENDENT POSITION FUNCTIONAL THEORY 

The evolution of electronic systems subjected to time- dependent potentials is in general 
too complicated to allow us to construct the functionals v[n;r, t], I^tM) etc., even in the 
simplest non-trivial case of a two-electron system. However, a simpler "position-functional 
theory" can be easily formulated, which is conceptually equivalent to the full-fledged theory 
and allows us to demonstrate explicitly all the main points of the theory. 

Our model is based on a two-particle system in one-dimension, with a time- dependent 
hamiltonian of the form 

H F (t) = X - [pi + pi + xj + xj] + - x 2 f - F{t)x x , (Bl) 

where X\, x 2 are the position operators of the two particles pi and p 2 the canonical mo- 
mentum operators, and F(t) is a time-dependent force, which acts only on particle 1. The 
two particles are subjected to a parabolic potential well, and interact with each other with 
a harmonic force with "elastic contant" k > 0. The idea is that x\ plays the role of the 
density operator; its expectation value Xi(t) is the time-dependent density, —F(t) is the 
external potential. As in TDDFT, one can show that the time-dependent position xi(t) and 
the initial state of the system at t — uniquely determine the force f(t) that produces it; 
but in this case the functional f[xi,t] can be explicitly constructed. 

For defmiteness, we start from an initial state described by the wave function 

^ (x 1 ,x 2 ) = Ce- x2 e-^~ k * 2 / 4 (B2) 

where X = (x± + x 2 )/2, x = x\ — x 2 , and C = ^172 1 is the normalization constant. This 
is the ground-state of the hamiltonian for F = 0. The time evolution of this state under the 
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full time-dependent hamiltonian is 



^x u x 2 ,t) = c e ^) e -[x-x cW P e2i [x-x c (t)]x c (t) 



x e -VT+2k{x-x c (t)] 2 /4 e i{x-x c (t)]x c (t)/2] 



(B3) 



where X c (t) and x c (t) are the solutions of the classical equations of motion 



X c {t)+X c {t) = F{t)/2 

x c {t) + (1 + 2k)x c {t) = F(t) (B4) 

with initial conditions X c (0) = X c (0) = and x c (0) = x c (0) = 0. The phase factor <fi(t) is 
the classical action (including the zero-point energy): 

<l>{t) = J + ^+^ t + f L(t')dt> , (B5) 

where 



o 



• 9 o x 2 r fl + 2k\ 9 F 

L = X 2 C - X 2 C +FX c + ^--( —— x\ + -x 



4 V 4 ) c 2 



(B6) 



is the classical Lagrangian. 

The solution of the equations of motion is 

•t 



1 r 

X c (t) = - sm(t - t')F(t')dt' 

2 Jo 



*.(«> - / 3il ' IV ™-' ,) ^ (0^ • (B7) 



The solution of the quantum mechanical problem is obtained by substituting Eqs. (1B7j) 
into Eqs. ( 1B3I) . ( 1B5I) . and ( 1B6I) . It is immediately evident that x c (t) and X c (t) are the 
expectation values of the quantum mechanical center of mass operator X = (xi + X2)/2 and 
relative position x = X\ — x<i respectively. The expectation values of X\ and x<i are given by 

x c {t) 



<^(t)|xi|V(t)> =x x {t) = X c (t) + 



2 



^{t)\x 2 m)) = x 2 (t) = X c (t) - ^ . (B8) 
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Our task is now to express the external force and the wave function as functionals of 
x\{t) - the "density" of our model. To do this, we observe that x 2 (t) is related to x\(t) by 
the classical equation of motion 

x 2 {t) + x 2 (t) = -k[x 2 (t) (B9) 

with initial condition x 2 (0) = x 2 (0) = 0. The solution of this equation, for given xi(t), is 

x 2 (t) = k — 1 v ^x x {t')dt' . (BIO) 



Jo VI + k 

From this we can express both X c (t) and x c (t) as functionals of Xi(t), and hence the whole 
time-dependent wave function ip{ x ii x 2, t) as a functional of Xi(t). Furthermore, the force /, 
which produces the evolution x\{t) is given by f(t) = x\{t) + x\{t) + k[xi(t) — x 2 (t)). Upon 
substituting the functional dependence of x 2 (t) on xi(t') in the expression for f(t) we obtain 
the force functional 

f[ Xl ;t} = Xl (t) + (1 + k) Xl (t) 



- k 2 f^^J^-^x l{ t>)dt>. 
Vl + k 



(Bll) 



Observe how the force is uniquely and causally determined by x\(t). Knowing the force we 
can construct the phase <f>(t) (Eq. (1B5I) ) as a functional of x\, by substituting F = f[xi] in 
the Lagrangian (1B6j) . Finally, we construct the internal action functional 

4)N = / (V'N;*]|^ t -#o|V'^i,*])^ 







I flxut'^ityt' , (B12) 

>/ 



where /[xi, t] is given by Eq. ( 1B11I) . 

We are now in a position to demonstrate explicitly the connection between the force and 
the functional derivative of the action. Namely, we can prove that 

5 -Prv) ■ (B13) 
Sx^t) I v ; 

where the state is described by the wave function (1B3[) . evaluated at time T and 

expressed as a functional of X\\t\ (the negative sign on the left hand side comes from the 



(5Ao[ii] ■ / / r i 
-/[^*] = -^-*\^r[x 1 ] 
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fact that the force enters the hamiltonian Hp (Eq. IB1I) with a sign opposite to that of the 
potential in Eq. (j3J.) 

The calculation is greatly simplified by the following two observations: (i) All the terms 
that involve an expectation value of X — X c or x — x c are obviously zero, and (ii) because 
x c (t) and X c (t) are solutions of the classical equation of motion, the variation of the phase 
4>(T) comes only from the variation of the force (regarded as a functional of of xi), and from 
the variation of X2(t) at the upper limit of integration. In this way we easily arrive at 

Jo 



- i ( i/j t [xi\ 



Sxi(t) 



x^dt' 



5xi(t) 

from which Eq. flB13j) follows at once. 

Similarly, we can show that the "Berry curvature" 2Sm ' 



(B14) 



Sx!(t') 



fapT l x l 



6x 1 {t) 



5x 2 {T) 5x 2 {T) 5x 2 {T) 5x 2 (T) 
SxiiV) Sx^t) 5xi(t) 5xi(t')' 

so making use of Eq. fIBlOl) we obtain 



is given by 

(B15) 



2Sm 



S^t\xi\ 
5xi{t') 



6il> T [xi] 



Sxi(t) 



,sm[Vl + k{t-f) 



(B16) 



VTTk 

from which the arbitrary time T has disappeared! Armed with this result it is an easy 
matter to verify that 



Sf[x i; t] 



XSx^t') 



5ip T [x 1 



5xi(t) 



(B17) 



<Jxi (f) 

for t > t', in agreement with Eq. (1251) . We can also verify the presence of the singularity at 
t = t' discussed in section [Til after Eq. (127]). Indeed, the functional derivative of the first two 
terms in the expression of our force functional, Eq. (IBllj) . gives 

^4 =-5(t-t') + (l + k)5(t-t'). (B18) 

Sx l(t) sing 

Notice that there is no term proportional to 5 in this simple model. 

Finally, we observe that the analogue of the xc potential - an xc force in this case - is 
F xc [xi] t] = f s [xi,t] — f[x\,t], where the non-interacting force functional f s [xi,t] is obtained 
from Eq. (IBlip simply by putting k = 0, so that 



fa[x u t] = Xi(t) + Xt(t) 



(B19) 
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and 

, / n , 2 f* sinfVl + Ht - t')] , .. 
F xc [xi,t} = -k Xl (t) + k 2 / [V ■_ L. ^ari(0 

Jo v 1 + « 

= -x 2 [si;t]}. (B20) 

Of course, F xc is nothing but the force exerted by the second particle on the first, expressed 
as a functional of the basic variable Xi(t). It is worth noting that the singular term Xi(t) 
has cancelled out in F xc . The singular part of the functional derivative of F xc is simply a 
5-function 

SF xc [xi,t] 



Sx^t') 



-k5{t-t'), (B21) 



sing 

in agreement with the discussion following Eq. (|28l) . 

Let us now demonstrate explicitly how, given the knowledge of the exact functional 
F xc [xi, t] (Eq. (lB20p ). one can calculate the evolution of Xi(t) within a "Kohn-Sham scheme" . 
First of all, we introduce the "Kohn-Sham hamiltonian" 

H,(t) = \ [Pi +Pl + x\ + xj] - F{t)x! - F xc [ Xl - t} Xl (t) , (B22) 

which describes two non-interacting particles, with an effective force F + F xc acting only on 
particle "1". Then we solve the time-dependent Schrodinger equation [id t —H s (t)]\ip s (t)) = 0, 
starting with the non-interacting ground-state 

Mx u x 2 ) = ^=e^ 2 e^ 2 , (B23) 

which clearly has the same expectation values of x\ and p\ as its interacting counterpart (1B2I) . 
The solution is the time- dependent "Kohn-Sham wave function" 

ijj s (x u x 2 ,t) = -^e-^i-^icWlVs^^i-xxeWJiicCt) 



(B24) 



x e x 2 /2 



where X\ c (t) is the solution of the equation of motion 

xUt) + xx(t) = F(t) + F xc [ Xl ; t] (B25) 
with initial conditions x lc (0) = Xi c (0) = 0, and 

Mt) = ~l + f L x {t')dt' , (B26) 
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with 



• 2 9 



M*) = -y - -y + ( F (t) + ^ c [xi c ; t])x lc (t) . 



(B27) 



Finally, we make use of the Kohn-Sham wave function to calculate the expectation value 
of X\\ 

tys{t)\x 1 \$ B {t))=x u {t) ■ (B28) 

The equation of motion (113250 for X\ c (t) can be rewritten explicitly as an integro- 
differential equation 



x lc {t) + (1 + k)x lc {t) = F(t) 



+k 2 



sin [VI + k(t - f )] 



x lc (t') 



(B29) 



with initial conditions Xi c (0) = Xi c (0) = 0. This equation can be solved by Laplace trans- 
formation. Denoting by Xi c (s) the Laplace transform of X\ c (t) we get 

k 2 



{s 2 + 1 + k)x lc (s) = F(s) + 



s 2 + 1 + k 



Xlc(s), 



(B30) 



and finally 



x lc (s) 



s 2 + 1 + k 



^ 2 + l)(s 2 + l + 2A;) 
1 1 



F{s) 



_s 2 + l s 2 + l + 2k 
Going back to the time domain we finally obtain 



F(s) 



(B31) 



i r* 

x lc (t) = - / sm(t - t')F(t ') 



+ 



sin[Vl + 2k{t - t')} 



F(t f ) 



(B32) 



2 Jo VTT2k 

which of course agrees with the exact solution Xx(t) = X c (t) +x c (t)/2 obtained from 
Eqs. (1B71). 

In TDDFT we do not have the luxury of knowing the exact xc force functional. But 
this example shows that, if we knew it, we could use it to predict the exact evolution of 
the density. So there is every reason to believe that a good approximation to the exact xc 
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potential, obtained by whatever means, would enable us to make good predictions for the 
time evolution of the density. 
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